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Abstract 

The dynamics of opinion formation in large groups of people is a complex non-linear phenomenon 
whose investigation is just at the beginning. Both collective behaviour and personal view play 
an important role in this mechanism. In the present work we mimic the dynamics of opinion 
formation of a group of agents, represented by two states ±1, as a stochastic response of each 
agent to the opinion of his/her neighbours in the social network and to feedback from the average 
opinion of the whole. In the light of recent studies, a scale-free Barabasi- Albert network has been 
selected to simulate the topology of the interactions. A turbulent-like dynamics, characterized by 
an intermittent behaviour, is observed for a certain range of the model parameters. The problem of 
uncertainty in decision taking is also addressed both from a topological point of view, using random 
and targeted removal of agents from the network, and by implementing a three state model, where 
the third state, zero, is related to the information available to each agent. Finally, the results of 
the model are tested against the best known network of social interactions: the stock market. A 
time series of daily closures of the Dow Jones index has been used as an indicator of the possible 
applicability of our model in the financial context. Good qualitative agreement is found. 
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I. INTRODUCTION 



Systems composed of many parts that interact with each other in a non-trivial way are 
often referred to as complex systems. The social relations between individuals can perhaps 
be included in this category. An intriguing issue concerns the role played by the topological 
structure of the social network in governing the dynamical behaviour of the system. 

Recent studies of the topological properties of interactions in different biological, social 
and technological systems has made it possible to shed some light on the basic principles 
of structural self-organization. A few examples include the food webs power grids and 
neural networks |2L |3J, cellular networks 4 1, sexual contacts 

World Wide Web 0, actor collaborations 2,0,0,0], the citation networks of scientists [l(3] 



and the stock market 11]. Although different in the underlying interaction dynamics or 
micro-physics, all these networks have shown a tendency to self-organize in structures that 
share common features. In particular, the number of connections, k, for each element, or 
node, of the network follow a power law distribution, P(k) ~ k~ a . Networks that fulfill this 
property are referred to as scale-free (SF) networks. In addition many of these networks are 
characterized by a high clustering coefficient, C, in comparison with random graphs [12]. 
The clustering coefficient, C, is computed as the average of local clustering, C*, for the ith 
node, defined as 

where Zi is the total number of nodes linked to the site % and yi is the total number of links 
between those nodes. As a consequence both Cj and C lie in the interval [0,1]. The high 
level of clustering found supports the idea that a herding phenomenon is a common feature 
in social and biological communities. 

Numerical studies on SF networks have demonstrated how the topology plays a funda- 
mental role in infection spreading [13] and tolerance against random and preferential node 
removal Q[. A detailed description of the progress in this emerging field of statistical 
mechanics can be found in the recent reviews of Refs. Q, Q. In the present work we 
investigate the implication of a scale-free topology in a stochastic o pini on formation model. 
Similar versions of this model have been tested in regular lattices jl tI . Q] and percolation 
clusters [19]. These models adopt a mean field approach where the interactions are extended 



between all the individuals in the lattice or cluster respectively. In contrast, our simulation 
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focuses on the role of short-range first-neighbour interactions for cases where the topological 
structure of the interactions is not trivial. 

In the next section we describe the model used for the simulation. In Sec. 11111 we show 
the results obtained numerically while in Sec. II VI we investigate the importance of failures 
in the network during the process of opinion formation while Sec. |V] the two state model 
is extended to three states and the numerical results are compared. In Sec. EU and I VI II we 
test the results of our simulations against the best known social network: the stock market. 
In particular the time series of average opinion fluctuations obtained with the model is 
compared with the time series of price variations for the Dow Jones index from 13/1/1930 
to 13/4/2004. The final section presents further discussion and conclusions. 



II. THE MODEL 



In the present work we investigate the opinion formation process of a group of N indi- 
viduals, represented by nodes on a SF network. The mechanism of opinion formation is 
simulated using stochastic heat bath dynamics with feedback. The opinion of each agent is 
of a Boolean type. That is, at each discrete time step, t, the opinion is represented by one of 
two possible states (or spin orientations), namely <Ji{t) = ±1, for the ith agent. A practical 
example could be the decision to buy, (Tj(t) = +1, or sell, (Tj(t) = —1, a stock in a virtual 
stock market. 

In order to mimic the scale-free network topology we make use of the Barabdsi- Albert 
model 0| . This is based on two main assumption: (i) linear growth and (ii) preferential 
attachment. In practice the network is initialized with mo disconnected nodes. At each step 
a new node with m edges is added to the pre-existing network. The probability that an edge 
of the new node is linked with the ith node is expressed by Tl(ki) = k^j J2j kj. The iteration 
of this preferential growing process yields a scale free network, P(k) ~ k~ a where a = 3. 

It is worth noting that the Barabdsi-Albert model cannot reproduce a high clustering 

coefficient. In fact, the value of this coefficient depends on the total number of nodes in the 
I I 

network |l5( and in the thermodynamic limit, iV — > oo, C — > 0. In principle the observed 
local clustering can play an important role in the opinion formation of groups of people, 



independent of their total number. In order to account for this, we introduce a 
in the growth process, namely the triad formation proposed by Holme and Kim 



urther step 



2jj. In this 
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FIG. 1: (Color online). Example of a scale- free network. The number of nodes is 500 with 
clustering probability 9 = 0.9 and mo = m = 2, so that each new node is linked twice. The number 
of nodes has been kept small in order to preserve the clarity of the plot. Note that, for such small 
networks, a large scale invariant range is obtained only if one considers the ensemble average over 
several realizations. This plot has been realized with the Pajek software [21 1 . 

case, if the new added node is linked with an older node, i, having other links, then with 
a certain probability, 9, the next link of the new node, if any remain, will be added to a 
randomly selected neighbour of node i. This method of introducing friends to friends, while 
preserving the scale-free nature of the networks, generates high clustering coefficients that 
do not depend on the number of nodes in the network. The only tunable parameter that 
changes the value of the clustering coefficient is the clustering probability 9. An example of 
SF network generated with this algorithm is shown in Fig. ^for 500 nodes *. 

Once the scale-free network has been built, we randomly assign the spin values, ±1, to 
every node and start the simulation of opinion formation. We neglect, in the first approxima- 
tion, the network dynamics. This is equivalent to assuming that the time scale for evolving 
the network is much longer that the time needed for people to make a decision. 

The dynamics of the spins follows a stochastic process that mimics the human uncer- 

* Another model for acquaintance networks, showing properties similar to the one presented in this work, 
has been proposed by Davidsen et al. [2^ . 
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tainty in decision making [18|, |l9j]. Values are updated synchronously according to a local 
probabilistic rule: <7j(t + 1) = +1 with probability pi and <7i(t + 1) = — 1 with probability 
1 — Pi. The probability Pi is determined, by analogy with heat bath dynamics with formal 
temperature k^T = 1, by 



*® = i + e-mty ( 2 ) 



where the local field, hit), is 



Ii(t) = a£(t)Ni 1 £ aj{t) + /^(t)r(t). (3) 

3=1 

The first term on the right-hand side of Eq. (jSJ) represents the time dependent interaction 
strengths between the node i and his/her Nj information sources, which are the first neigh- 
bours in the network. The second term instead reflects the personal reaction to the system 
feedback, that is the average opinion, 

1 N 

K*) = T7£^(*), (4) 
JV 3=1 

resulting from the previous time step. The terms and f]i(t) are random variables uni- 
formly distributed in the interval (-1,1) with no correlation in time nor in the network. They 
represent the conviction, at time t, with which agent i responds to his/her group (common 
for all the agents) and the global opinion of the network respectively. The strength term, 
a, is constant and common for the whole network, while hi is specifically chosen for every 
individual from a uniform distribution in (0,k) and are both constant in the dynamics of the 
system. By varying the parameter k we can give more or less weight to the role of feedback 
in the model. The strength coefficients a and hi in the local field, ij, characterizing the 
attributes of the agents, play a key role in the dynamics of the model. They represent the 
relative importance that each agent of the network gives, respectively, to his/her group and 
to the variation of the average opinion itself. 

III. NUMERICAL SIMULATIONS 

At first we investigate the importance of the group strength a by fixing k — a. In 
this case the dynamical behaviour is similar to that found in the stock market context in 
Refs. For a < 1 the resulting time series of average opinion is largely uncorrelated 

Gaussian noise with no particularly interesting features, as illustrated in Fig. |2(i). 
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FIG. 2: Time series of the average opinion, r, for different values of the group interaction strength 
parameter a: (i) a = 0.8, (ii) a = 1.5, (iii) a = 1.8 and (iv) a = 2.3. The parameters used for the 
simulations are: N = 10 4 nodes, clustering probability 6 = 0.9, initial nodes and links per new 
node mo = to = 5 and we take the upper bound of the distribution of personal response strengths 
equal to the group interaction strength, that is k = a. The results involve 10 realizations of the 
scale free network each displayed for 5000 time steps. For values of a greater than 1 a turbulent-like 
state, characterized by large fluctuations, starts to appear in the process of opinion formation. The 
clustering probability = 0.9, related to the triad formation in the network, fixes the clustering 



coefficient to C « 0.39. This value is similar to that found for many real systems 
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As soon as we exceed the value of a ~ 1 a turbulent-like regime sets in, characterized by 
large intermittent fluctuations, as illustrated in Fig. |2[ii — > iv). These large fluctuations, 
or coherent events, can be interpreted in terms of a multiplicative stochastic process with a 



weak additive noise background 



18, 



231 ] . For a > 2.7 we observe that the bursts of the time 



series begin to saturate the bounds — 1 < r < 1. 

In Fig. El we plot the probability distribution functions (PDFs) Associated with the time 
series of Fig. |21 The large fluctuations, for a greater than pa 1, are reflected in the fat tails 
of the relative PDFs. Decreasing the value of a, and so the number of coherent events, the 
PDF converges to a Gaussian distribution generated by a random Poisson process. 

The personal response to the change in the average opinion also plays an important role 
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FIG. 3: PDFs of the time series relative to Fig. |2J The shapes of the distributions converge to a 
Gaussian for small values of the group interaction strength a = k. A Gaussian distribution is also 
plotted for comparison. All the PDFs in this paper are obtained over 50 realizations of the SF 
network. In order to compare the fluctuations at different scales, the time series in the plot have 
been normalized according to r(t) — ► r ^ r , where r and a denote the average and the standard 
deviation over the period considered respectively. 

in the turbulent-like regime of the simulation. In order to study the impact of this term 
on the dynamics we change the parameter k while keeping a fixed at 1.8. The results are 
summarized by the PDF plots in Fig. 0] For k = the behaviour of the time series is still 
turbulent- like, underlying how the network group interaction is, in reality, the only crucial 
factor for the appearance of coherent events. As expected, incrementing the value of k leads 
to a progressive crossover toward a noise regime. It is important to notice how this regime is 
reached for k > 10a. The group interactions continue to play an important role even when 
the average value of hi is large compared to a. 

In order to test the relevance of the network structure on the process of opinion formation, 
the previous simulations have been repeated, with a large number of nodes, N, and k — a, 
for different values of the clustering parameter, 9, and the node-edge parameter, m. While 
varying 6, does not lead to any substantial difference in the dynamics of the model, the 
increase of the average number of links per node, k = 2m, has a dramatic effect in the 
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FIG. 4: The importance of the personal response, related to the global opinion strength parameter 
K is shown by the change of shapes of the PDFs for group interaction strength a = 1.8. For large 
values of k the time series of global opinion approaches Gaussian noise. The time series of r has 
been normalized - see the caption of Fig. |3 

turbulent-like phase, as shown in Fig. El Here the kurtosis, (K r = (r 4 )/(r 2 ) 2 , where (...) 
denotes the temporal average), of the time series of the average opinion, used to quantify 
the deviation from a noise regime, is plotted against m. It is evident that an increase in 
the average number of links per node gives rise to more turbulence characterized by larger 
fluctuations and broader tails in the PDF. Large scale synchronizations are more likely to 
occur for large m. This behaviour is intrinsically related to the model of Eqs. (0) and (j3J). In 
fact, the turbulent-like regime is a consequence of the random fluctuations of the interaction 
strengths between agents around a bifurcation value separating the ordered and disordered 
phase. 

If we take the thermodynamic limit, where iV — > oo and m — > oo , then the coupling 
strengths between agents can be approximated well by the average strength over all the 



network and a mean fie 
model. Krawiecki et al. 



d approach becomes appropriate to describe the dynamics of the 
proposed the following map 



r{t + l)=A£{t)r(t) + hTi{t), (5) 
as a mean field approximation of a stochastic dynamical system similar to the one used in 
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FIG. 5: Dependence of the kurtosis, defined as K r = (r 4 )/(r 2 ) 2 , where (. . .) denotes the temporal 
average, as a function of the node-edge parameter m. For a Gaussian noise process K r = 3 while 
for K r > 3 large deviations from the average start to appear. The final value for each m has 
been obtain after the average over 50 configurations of the network. The calculations show an 
exponential increase for K r . 

the present work. Here A and h are coupling coefficients and and r](t) random numbers 



in the interval (-1,1). The map of Eq. © is a generic mode 



24, 



3, 



for on-off intermittency and 
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attractor bubbling extensively studied in chaos theory 

It is also worth pointing out that an increase of k is related to a decrease in the average 
path length between nodes; that is, the network "shrinks" and becomes more compact. In 
relation to our previous discussion, the more compact the network is the more the dynamics 
of our system approaches to the mean field approximation. It becomes easier for the agents 
to synchronize. This characteristic of compactness, referred to as the small world effect 
Q Q| , is actually very common in both real and artificial networks. 

We further investigate the importance of the SF network topology and the the small world 



effect in our model by performing a numerical simulation of the same system but using a 
random network (RN) or random graph as the underlying topology. Given a fixed a 
nodes, N, a RN is defined by the probability p that two nodes are linked together [ 
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FIG. 6: Left: comparison between the PDFs of our model obtained on a SF network and on a 
RN with number of nodes N = 10 4 and average links per node k = 10. For the SF network the 
parameters used are m = m,Q = 5 for the links of each new node and 9 = 0.9 for the clustering 
probability while for the RN p = 10/ N. From a statistical point of view the characteristic features 
of the PDFs have their origin in the model dynamics as opposed to the fine features of the network. 
Right: Dependence of the opinion fluctuations on the parameter p on a RN. The parameters used 
for the dynamics are a = 1.8 and n = a for the group and global opinion response respectively. 

In this case k = pN and, moreover, there exists a critical value, p = p c ~ 1/N, for which the 
the network undergoes a topological phase transition where it moves from a phase where it 
is composed of a collection of small, disjoint, sub-networks to a phase where a giant cluster 
emerges b Random networks, while preserving small world properties, have a Poisson degree 
distribution [15], P(k), and small clustering coefficients. As previously mentioned, we make 
use of the RN to test the robustness of our model with respect to the topology used and 
to learn about the most important properties relevant to the dynamics. In order to do so, 
we fix the number of agents and the average number of links for the SF and RN, namely 
N = 10 4 and k = 10. Then we perform independent numerical simulations on the two 
topologies. Note that for the RN, k = 10 requires p = 10/iV, that is ten times greater than 

t Note the analogy between the random network theory and the standard percolation theory on a lattice 
where the structural properties of the system are studied as a function a percolation probability. 
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the percolation threshold. The results, shown in Fig. |H1 (left) demonstrate how the dynamics 
of the two systems are largely equivalent under the adopted constrains. In Fig. El (right) we 
also show the dependence of the dynamics on the parameter p for the RN. At the critical 
threshold, that is the value of p for which a giant cluster appears, there is still no trace of 
turbulent- like activity giving rise to fat tails. Yet, in this case each agent has, on average, 
just one link and there cannot be any small world properties. 

These results confirm that the critical topological characteristic leading to herding be- 
haviour in the framework of stochastic opinion formation is the presence of mean field effects 
enhanced by small-world structure. The more information (links) that an agent has, the more 
likely it is for him/her to have an opinion in accord with other agents. 

In the next section we extend our model in order to include indecision in the process of 
opinion formation. 

IV. THE INFLUENCE OF INDECISION 

We now extend our model in order to include the concept of indecision. In practice a 
certain agent i, at a time step t, may take neither of the two possible decisions, o~i = ±1, but 
remain in a neutral state. Keeping faith to the spirit of the model, we address this problem 
introducing an indecision probability, e: that is the probability to find, at each time step, 
a certain agent undecided. This is equivalent to introducing time dependent failures in the 
structure of the network by setting a = 0. 

Focusing on the turbulent-like regime, the shape of the PDF in the opinion fluctuations 
changes according to different concentrations of undecided persons. The results of the sim- 
ulations, in Fig. [7[ show how the dynamics of the model move from an intermittent state 
for e = toward a noise state for e ~ 0.6. The convergence to a Gaussian distribution is ob- 
tained only for quite high concentrations of undecided agents at about 60%. The robustness 
of the turbulent-like behaviour is related to the intrinsic robustness of SF networks against 
raudom failu.es Q. In fact, became there is a la rg e absolute uuruber of poo rly couueeted 
nodes, related to the power law shape of P(k), the probability of setting one of them to 
inactive is much higher compared to the "hubs" that are relatively rare. 

We can claim that, in large social networks governed by stochastic reactions in their 
elements, large fluctuations in the average opinion can appear even in the case in which a 
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FIG. 7: Transition from coherent bahaviour, indecision probability e = 0, to noise using a random 
selection for the inactive agents. For e ~ 0.6 we reach a noise-like behaviour. The parameters used 
in the simulation are: N = 10 4 nodes, 9 = 0.9 for the clustering probability, m = ttiq = 5 for the 
links of each new node, a = 1.8 and n = a for the group and global opinion response respectively. 

large part of the network is actually "inactive" provided that the structure is scale free and 
the indecision is randomly distributed. The existence of large hubs provides for the survival 
of extended sub-networks in which synchronization can give rise to coherent events. The 
structure of the network itself supplies the random indecision. 

Now we address the question of how the dynamics may change if we do not choose 
randomly the inactive nodes but we target the nodes having the most links. What we 
do in practice is to sort the nodes according to their number of links and then deactivate 
the nodes having the largest number of links in decreasing order. Fig. |S] illustrates how 
the fragmentation process is much faster and the noise regime is reached already when 
only the 10% of the hubs are deactivated. As emphasized in Ref. the hubs have a great 
importance in the structural properties of SF networks and specifically targeting these nodes 
can lead to sudden isolation of a large fraction of the nodes of the network. 
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FIG. 8: In this simulation we progressively turn off the largest hubs in the network. Once we have 
turned off about the 10% of agents, N = 10 4 , the coherence in opinion formation disappears. The 
parameters used in the simulations are the same as in Fig. [7| 

V. AGENT INDUCED INDECISION: THE THREE STATE MODEL 

In the previous section we introduced random and targeted failures in order to study the 
response of the system to changes in the network topology. In a real social network the 
reason behind the indecision of a person follows much more complex rules and can depend 
on different factors as, for example, unsatisfactory information obtained by his/her sources. 
As seen from Eq. (jBJ), the opinion of each agent depends on the poll of his/her network links. 
Suppose now that the agent % has iVj neighbours where iVj/2 of these share the opinion +1 
while the remaining iVj/2 share the opposite opinion. In this case, unless we give specific 
weights to each node, the agent i will not have an easy task in choosing one of the two 
possible positions because of a lack of popular consensus. Based on this idea derived from 
common sense, we can extend our two state model by introducing an induced indecision 
'probability, fi, dependent on the information available to the agents at each time step. In 
particular we define the global opinion of the neighbours of the ith node as Sj(t) = J2f=i 
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and the indecision probability for the ith node at time t 



(6) 



where the indecision probability width, <j, is a parameter of the model and q a normalization 
constant that depends just on the structure of the network. It calculated at the beginning of 
the simulation by imposing Y^^i_ft. fais, 0) = 1, i.e. the sum of the indecision probabilities 
over all possible global opinions to be one. The model of Eq. (jBJ) assumes a Gaussian 
probability, centered in s, = 0, for the distribution of indecision of the ith agent. That is, 
the probability of having this agent in a state with Oi = is greater when there is not a 
large agreement in the opinion of the his/her sources. 

The analysis of the time series generated by the three state model does not present any 
relevant difference if compared with the two state model with the same parameters, Fig. 

We also plot the PDF for the number of inactive agents, N s (t), during the simulation, 
Fig. E3 It is interesting to notice how this distribution is not Gaussian distributed around 
the average but it is skewed on one side. Moreover, only a small fraction of agents is 
undecided, of the order of 10/15 %. This is consistent with the observation that in opinion 
polls most of the participants actually indicate an opinion. 



VI. POSSIBLE APPLICATION: OPINION FORMATION AND THE STOCK 
MARKET 



The model for opinion formation discussed thus far can be tested against the best known 
real social network: the stock market. The main idea is to compare our results with some 
stylized facts concerning the price time series, P(t) and, in particular, with the properties 
of the logarithm of the price fluctuations, or returns, R(t) = In P(t + 1) — lnP(t). In fact 
some characteristic features are independent of the particular market and can be considered 
as universal y. Moreover the returns show an intermittent behaviour, reminiscent of 



hydrodynamic turbulence (30 



31 



3, 



331 ] . also characterized by power law tails in the PDF. 



In this case the large coherent events are related to crashes or other anomalous variations 
of price. 

If we assume that the variation of price is directly proportional to changes in demand 
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FIG. 9: (a) A window of the normalized time series generated by the two-state model with 
parameters N = 10 4 nodes, 9 = 0.9 for the clustering probability, m = rriQ = 5 for the links of each 
new node, a = 1.8 and n = a for the group and global opinion response respectively, (b) Window 
of the normalized time series generated by the three states model with the same parameters as (a) 
and indecision probability width <r = 1. (c) Comparison between the PDFs generated by the two 
and three-state models with the aforementioned parameters obtained over 50 realizations of the 
SF network. No relevant differences can be observed. 



and supply, 



dP 



(7) 



where c p is proportional to the average opinion, r(t), then the returns are proportional to 
the average opinion R(t) ~ r(t). Using this assumption, we compare the time series of 
average opinion generated by the two state model against the time series of daily closures of 
the Dow Jones index. The data set spans the range 13/1/1930 to 13/4/2004 for a total of 
18645 samples. In Fig. El a comparison between the two PDFs is shown. The similarities 
between the model and the Dow Jones is remarkable. Both distributions have a leptokurtic 
shape and, in particular, they are described by power law tails, expressing the turbulent- like 
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FIG. 10: PDF of the number of inactive agents, cr^(t) = 0, during the simulation of the three state 
model. The parameters used are the same as used in Fig. |HI 



dynamics of the time series*. Note that, in contrast to the self-organized model for stock 
market dynamics proposed by Bak et al. 24|, here the price feedback is not an essential 
ingredient for the reproduction of the correct shape of the distribution. Rather it is the 
herding behaviour that plays the main role, as observed from Fig. HJ 

The similarities between the artificial time series generated by the virtual social network 
and the stock market extend beyond the fat tails in of PDF of the fluctuations to temporal 
correlations. It is well known that the stock market returns have negligible correlations on 
daily intervals while the volatility, u, denned as their absolute value, have a slow power law 
decrease as a function of the time lag. This phenomenon is known as volatility clustering . 
In order to make a comparison with our model we make use of the autocorrelation function, 
p. For a time series of L samples, Xi for i = 1, L, this is defined as 



P(r) 



j=i \ x j 



x){x 



3+t 



X) 



J2j=i( x j 



x)< 



(8) 



* The problem of the actual shape of the PDF for the stock market returns is still a matter of debate in 
the economics community HQ A solution to this problem would be of a great interest, especially 
for the practical application of option pricing. 
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FIG. 11: Comparison between the PDF of our model and the time series of the Dow Jones index 
from 13/1/1930 to 13/4/2004. The parameters of the model used to reproduce the PDF in the 
plot are: N = 10 4 nodes, 9 = 0.9 for the clustering probability, m = too = 5 for links of each new 
node, a = 1.8 and k = a for the group and global opinion response respectively. A Gaussian is also 
superimposed in order to emphasize the fat tails. 

where r is a time delay and x represents the average over the period under consideration. 
The autocorrelation has been computed both for the returns and for the volatility. While the 
time series of returns generated by the model and the Dow Jones index have an equivalent 
behaviour, Fig. (top), the same similarities do not hold for the volatility, Fig. 1121 fbottomb 
We observe a qualitatively different correlation: while for the market we observe a power law 
behaviour, the memory in the time series gen erated by the model decays exponentially like a 
short-range correlated random processes j3fj. This second point illustrates how non-trivial 
memory effects in the stock market cannot be taken into account by a simple heath bath 
dynamics. 

In Fig. (bottom) we also reproduce the autocorrelation function for the model pre- 
sented in Ref. |19.j. In this model an heat bath dynamics, similar to the one used in the 
present simulations, is applied to dynamical percolation clusters, used as a paradigm for 
agents aggregation. The temporal evolution of the clusters, which size follows a power law 
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FIG. 12: Autocorrelation functions for the fluctuations r (top) and the volatility v (bottom). 
The parameters used to produce the analyzed set are: N = 10 4 nodes, 9 = 0.9 for the clustering 
probability, m = m,Q = 5 for the links of each new node, a = 1.8 and k = a for the group and 
global opinion response respectively. 

distribution, is related to a forest-fire dynamics in which some potential traders are attracted 
in the market by other already active traders while, at the same time, some of them may 
temporally quit the trading. Large fluctuations in the price changes are due to the synchro- 
nization in the of the larger clusters in the market at a particular time. The main qualitative 
difference between this model and the one presented so far is that the former presents a de- 
cay rate much closer to that of the real market. At this point it is important to underline 

ated to the network dynamics. While 



that the main difference between the two models is re 
in the present simulation the network is fixed, in Ref. 19| the interaction between agents are 
time dependent and localized in separate clusters. We can argue that the dynamics of the 
networks and, in particular, the clustering of agents in different sub-networks can play an 
important role in the correlation properties of the stock market volatility. In reality, this fact 
appears quite natural if we use the autocorrelation function, defined in Eq. ©, in order to 
estimate the degree of memory in a process. If, for example, the variable under investigation 
is the sum over many independent Markovian processes, as in Ref. |l9(, then the resulting 
autocorrelation is given the convolution of the common exponential decay, oc e _/3r , with the 
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distribution of the decaying rates, g(fl), 

p(r) oc / g(P)e~^dp. (9) 



According to the shape of this distribution, the observed macroscopic variable can show a 
behaviour characteristic of a long memory processes, like the 1/f Fourier spectrum (3^. 
Power law tails in the probability distribution function, p(r) = r~ 7 , are produced from the 
distribution g(f3) = T( / ~f)^ 1 (3' y ^ 1 , where T is the gamma function and 7 a generic real expo- 
nent This fact strengthens the idea that the stock market is organized in a hierarchy of 
sub-networks where each of them can be considered, from a physical point of view, at local 
equilibrium. For time periods shorter than the typical time scale necessary for the networks 
to evolve, the only link between the sub-systems composing the market is the feedback com- 
ing from the price history. This idea is closely related to the concept of subordination used 
in probability theory |22j. The superposition of distributions, as a possible explanation of 



fat-tailed processes, has been proposed recently by Beck [29J in the context of hydrodynamic 
turbulence and then extended also to other systems j^j including the stock market 41 1. 



VII. MULTIFRACTAL ANALYSIS 



Financial time series present an inherent multifractality 42]. In the past few years the 



work of many authors 
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43, 



44. 



45 



has been addressed to the characterization of 
the multifractal properties of financial time series, and nowadays multifractality can be 
considered as a stylized fact. In order to study the multifractal properties of our model we 
use the generalized Hurst exponent |47(, H(q), derived via the q— order structure function, 



S q (T) = (\x(t + T)-x(t)\«) T <XT* B M, 



(10) 



where x(t) is a stochastic variable over a time interval T and the time delay, r. The 
generalized Hurst exponent, defined in Eq. QlOjl . is an extension of the Hurst exponent, 
H, introduced in the context of reservoir control on the Nile river dam project, around 



1907 142 



481 ]. This technique provides a sensitive method for revealing long-term correlations 



in random processes. If H(q) = H for every q, the process is said to be monofractal and 
H is equivalent to the original definition of the Hurst exponent. This is the case of simple 
Brownian motion or fractional Brownian motion. 
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FIG. 13: Structure function exponents for the Dow Jones index and our model. A deviation from 
a linear behaviour is evident. The hypothetical spectrum of a ID Brownian motion is also shown 
for comparison. 



If the spectrum of H(q) is not constant with q the process is said to be multifractal. From 
the definition it is easy to see that the function H(l) is related to the scaling properties 
of the volatility. By analogy with the classical Hurst analysis, a phenomenon is said to be 
persistent if H(l) > 1/2 and antipersistent if H(l) < 1/2. For uncorrelated increments, as 
in Brownian motion, H(l) = 1/2. In Fig.^Ja comparison is shown between the multifractal 
spectra of the model and the Dow Jones index obtained from the price time series. It is 
clear that both processes have a multifractal structure and the price fluctuations cannot be 

n 

associated with a simple random walk as in the classical efficient market hypothesis [49J . 



VIII. DISCUSSIONS AND CONCLUSIONS 



In the present work we have introduced a two state model of opinion in order to simulate 
the complex dynamics of opinion formation in a group of individuals. The decision updating 
is governed by a stochastic heat-bath dynamics that mimics the reaction of each person to 
his/her specific sources of information as governed by the network neighbours and to the 
average opinion of the whole group. Particular emphasis has been given to the topology of 
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the interactions between agents, where a Barabdsi- Albert scale-free network has been used 
to simulate the links between them. The choice of this particular network is motivated by a 
series of recent studies on social aggregation 1 15, 16] but, as we have shown in Sec. 11111 its 

\ nnnr 

use is not essential for the appearance of coherent events. As in other studies [17L I18L I19t| . we 
find a range in the parameter space in which the fluctuations of opinion have a non-trivial 
turbulent- like dynamics. The results of the simulations show that the most important factor 
determining the appearance of large fluctuations, is the synchronization of large parts of the 
network. As discussed in Sec. IIII[ this feature plays an important role even in the case in 
which the personal opinion is relatively strong. As a consequence large coherent events are 
more likely to occur when the average number of links per agent is larger. 

The topology of the interactions also plays a key part in the dynamics of the model. 
In fact, introducing inactive agents and spreading the undecided agents randomly on the 
network, does not spoil the turbulent-like state even for high concentrations of "gaps" , up to 
approximately 60% of agents. This is a consequence of the implicit robustness of SF networks 
against random failures. If instead of selecting randomly the undecided individuals we aim 
directly to the "hubs" of the network then the situation changes. In this case the network 
is disaggregate, composed of very small sub-networks and isolated nodes. Synchronization 
cannot significantly effect the resulting global opinion and the time series approximates 
Gaussian noise. We also introduce, in Sec. |V] a three state model. While the dynamics 
does not significantly differ from the two state model, we find a persistence of opinion with 
a sharp upper limit in the number of undecided agents. In Sec. EH we test the results of 
the simulations against a time series of daily closures for the Dow Jones index. The stock 
market, in fact, can be considered as the most studied network of social interactions. The 
results show a very good agreement with some stylized facts of the financial market like the 
broad tails in the PDFs, temporal correlations and a multifractal spectrum. We also notice 
an interesting discrepancy in the autocorrelation function for the volatility. Comparing the 
present results with those obtained in Ref. we conjecture that the persistence in the 
volatility memory can be explained by considering the market as constituted by sub-systems 
at local equilibrium and weakly interacting with each other. It will be interesting to explore 
this conjecture in a quantitative manner in a further investigation. 
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